THE INTEGRALS IN GRADSHTEYN AND RYZHIK. PART 4: 
THE GAMMA FUNCTION. 



VICTOR H. MOLL 



Abstract. We present a systematic derivation of some definite integrals in 
the classical table of Gradshteyn and Ryzhik that can be reduced to the gamma 
function. 



1. Introduction 

The table of integrals [2] contains some evaluations that can be derived by ele- 
mentary means from the gamma function^ defined by 



The convergence of the integral in (|l.ip requires a > 0. The goal of this paper 
is to present some of these evaluations in a systematic manner. The techniques 
developed here will be employed in future publications. The reader will find in [1] 
analytic information about this important function. 

The gamma function represents the extension of factorials to real parameters. 
The value 



(1.1) 




(1.2) T{n) = (n- 1)!, for n e N 

is elementary. On the other hand, the special value 




is equivalent to the well-known normal integral 



(1.4) 




The reader will find in [T] proofs of Legendre's duplication formula 



(1.5) 




that produces for a; = m e N the values 



(1.6) 




This appears as 3.371 in [2\. 
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2. The introduction of a parameter 

The presence of a parameter in a definite integral provides great amount of 
flexibility. The change of variables a; = /i< in (jl.ip yields 



(2.1) r(a) / t°-ie~^*dt. 

Jo 

This appears as 3.381.4 in and the choice a = n + 1, with n e N, that reads 



(2.2) / fe-^"' dt = 111 ^1-''-^ 

Jo 

appears as 3.351.3. 

The special case a ~ m + ^, that appears as 3.371 in [2], yields 

(2.3) rr-ie--*dt = #i^^— i 
^ ' io 22™ ml ^ 

is consistent with (|1.6p . 

The combination 

(2.4) / ' d^^^^^r(i-p), 

that appears as 3.434.1 in [2\ can now be evaluated directly. The parameters are 
restricted by convergence: fi, > and p < 1. The integral 3.434.2 



e - e ly 
(2.5) / dx = ln-, 



X 



is obtained from (|2.4p by passing to the limit as p ^ 0. This is an example of 
Frullani integrals that will be discussed in a future publication. 

The reader will be able to check 3.478.1: 



/ a;''-iexp(-/ia;P) dx = -^""/f F 
7o P 



(2.6) 

and 3.478.2: 

(2.7) / a;''-Ml-exp(-//a;P)] dx = ---p-''/Pr - 
Jo \P\ \P 

by introducing appropriate parameter reduction. 

The parameters can be used to prove many of the classical identities for r(a). 
Proposition 2.1. The gamma function satisfies 

(2.8) r(a + l) =ar(a). 
Proof. Differentiate ()2.ip with respect to /i to produce 

/•oo /-oo 

(2.9) O^a/^'^^W t°^ie"^*di-//° / i^e"'^* dt. 

Now put p = 1 to obtain the result. □ 
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Differentiating (jl.ip with respect to the parameter a yields 

/•oo 

(2.10) r'(a) = / x^-^e-"" \iYxdx. 

Jq 

Further differentiation introduces higher powers of In a;: 

/■oo 

(2.11) r(")(a)=/ x^-^e-^ (Inx)" dx. 

Jo 

In particular, for a = 1, we obtain: 

/•OC 

(2.12) / (lnx)"e-^dx==r(")(l). 



The special case n — 1 yields 
(2.13) 



/•OO 

/ e"^ lna;da; = r'(l). 
Jo 

The reader will find in [I], page 176 an elementary proof that r'(l) = —7, where 

n ^ 

(2.14) 7 := lim V - - ln?i 



is Euler's constant. This is one of the fundamental numbers of Analysis. 
On the other hand, differentiating (|2.1p produces 

(2.15) j x'^-^e-^''' {\nxf dx= [^^"^(a)] , 

that appears as 4.358.5 in [5]. Using Leibnitz's differentiation formula we obtain 

(2.16) / x°-ie-^"(lna;)" fix = Ai"''5](-l)''(, I (ln/i)''r("-'=)(a). 
In the special case o = 1 we obtain 



(2.17) 



/•OO 1 '^'^ / \ 

/ e-^-(lnx)" dx = -y^{-lfr] (In^)" r("-'=)(l). 



The cases n = 1, 2, 3 appear as 4.331.1, 4.335.1 and 4.335.3 respectively. 

In order to obtain analytic expressions for the terms r'^"^(l), it is convenient to 
introduce the polygamma function 

(2.18) ip(x) ^ -^\nT{x). 

ax 

The derivatives of ip satisfy 

(2.19) ^/>(")(a;) = (-l)"+in! ({n + 1, x), 
where 

00 

n=0 ^ ^' 

is the Hurwitz zeta function. In particular this gives 
(2.21) ^l;^"\l) = {-l)"+^nl({n+l). 
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The values of r(")(l) can now be computed by recurrence via 

(2.22) r("+i)(i) = (f\r^''^i)i^^"-''\i), 

k=o ^ ^ 

obtained by differentiating r'(a;) = 'ip{x)T{x). 

Using (|2.f 9p the reader wiU be able to check the first few cases of (|2.15p . we 
employ the notation 6 — ip{a) — In/i: 

Jo 

poo yfn) 

/ a;''-ie-^^ln3a;da; = {6^ + 3C{2, a)6 - 2C{3, a)} , 

Jo M° 

f x^-h-f^'ln^xdx = ^-^{^* + 6C(2,a)(52-8C(3,a)(5 + 3C^(2,a) + 6C(4,a)} 
Jo 

These appear as 4.358.2, 4.358.3 and 4.358.4, respectively. 

3. Elementary changes of variables 

The use of appropriate changes of variables yields, from the basic definition (II. 
the evaluation of more complicated definite integrals. For example, let x = io 
obtain, with c = ab — I, 

1"°° 1 
(3.1) / t''eyi^{-t'')dt = -T 

Jo b 

The special case a — 1/6, that is c = 0, is 
(3.2) 



c+ 1 



that appears as 3.326.1 in [2]. The special case 6 = 2 is the normal integral (|1.4p . 
We can now introduce an extra parameter via t — s^^^x. This produces 

(3.3) x"exp(-sx'') dx = 

with m = ab—1. This formula appears (at least) three times in [5]: 3.326.2, 3.462.9 
and 3.478.1. Moreover, the case s = 1, c = (m + l/2)n — 1 and h = n appears as 
3.473: 

(3.4, a^^^ 

Jo 2'"n 
The form given here can be established using (ll.6|l . 

Differentiating p.3p with respect to the parameter m (keeping in mind that 
a = (to + l)/6), yields 

(3.5) / x^e""^' lna;da; = ^J^^[i/;(a)-lns]. 
In particular, if 6 = 1 we obtain 

(3.6) j x^e""^ \nxdx= ^^"^,^1^^ [7/>(to + 1) - In a] . 
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The case m = and b — 2 gives 

(3.7) / e""^' ln2;d2;==--^(7 + ln4s), 

where we have used V'll/S) = ^7 ~ 21n2. This appears as 4.333 in [2]. 
An interesting example is = m = 2. Using the values 

(3.8) r(f) =\/^/2and =2-21n2-7 
the expression (|3.5p yields 



(3.9) / x^e""^" Inxda;^ — (2-ln4s-7),/^. 

h 8s \ s 

The values of '0 at half-integers follow directly from (jl.Sp . Formula (|3.9p appears 
as 4.355.1 in [2]. Using (|3.5p it is easy to verify 



(3.10) 

and 

(3.11) 



io 2(2^^)" 



for n € N. These appear as, respectively, 4.355.3 and 4.355.4 in [2]. The term 
[2n— 1)!! is the semi-factorial defined by 

(3.12) (2n - 1)!! = (2n - l)(2n - 3) • • • 5 • 3 • 1. 
Finally, formula 4.369.1 in 2J 

(3.13) r x'^-^e-''" {i,{a) - Inx] dx = ^^"^ ^'^^ 



can be established by the methods developed here. The more ambitious reader will 
check that 

\ x^'\'^- {[In a; - \^{n)f - (n)) dx = {[In^. - l^lj{n)f + ^^P^n] 

Jo M 

that is 4.369.2 in [2]. 



We can also write p.Sp in the exponential scale to obtain 



(3.14) y_^te-expHe-) = (V^ (y j - In ^ 
The special case b = m = I produces 

(3.15) r te' exp (-se*) dt = _ll±M 

that appears as 3.481.1. The second special case, appearing as 3.481.2, is b 

2, m = 1, that yields 

(3.16) r te*exp(-.e2*) dt^-^^^^^. 
This uses the value V'(l/2) = -(7 + 21n2). 
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There are many other possible changes of variables that lead to interesting eval- 
uations. We conclude this section with one more: let x = e* to convert (jl.ip into 

/oo 
exp(-e^) e'^'^dx = r(a). 
-OO 

This is 3.328 in [2]. 

As usual one should not prejudge the difficulty of a problem: the example 3.471.3 
states that 

(3.18) ^ a;-^-i(a-x)^-^e-'^/^dx = /3-^a''-ir(/i) exp(-f) . 

This can be reduced to the basic formula for the gamma function. Indeed, the 
change of variables t = P/x produces 

(3.19) / = /3-^a'^-M {t- p/af-'^ e'Ut. 

Jp/a 

Now let y — t — (j/a to complete the evaluation. The table writes fj, instead 
of a: it seems to be a bad idea to have /i and u in the same formula, it leads to 
typographical errors that should be avoided. 

Another simple change of variables gives the evaluation of 3.324.2: 

/°° 1 / 1 

n \2n 

The symmetry yields 

(3.21) 1^2 £-(^-6/^)'" dx. 

Jo 

The change of variables t = h/x yields, using 6 > 0, 

(3.22) J = 2fe /"e-^*-"/*)"^. 
The average of these forms produces 

(3.23) / = / (l + 4) dx. 



Finally, the change of variables u = x — b/x gives the result. Indeed, let u ^ x — b/x 
and observe that u is increasing when & > 0. This restriction is missing in the table. 
Then we get 



(3.24) ^ = 2 / e"" du. 

This can now be evaluated via v — 



GAMMA FUNCTION 



7 



Note. In the case 6 < the change of variables u — x — b/ x has an inverse with 
two branches, sphtting at a; = \J~-h. Then we write 

/"OO 

(3.25) / 2 / e-(^-''/=^)'" 



da; + 2 / e-(^-^/^)'" dx. 

The change of variables u = a; — 6/a;is now used in each of the integrals to produce 

uexp(— u^") du 



(3.26) / = 2 



The change of variables z ~ \/u^ + Ab yields 

/•OO 

(3.27) ^ = 2 / exp(-(z2 -46)") . 

Jo 

We are unable to simplify it any further. 

4. The logarithmic scale 
Euler prefered the version 

(4.1) T{a) = J^ i^^u) 
We will write this as 

(4.2) r(a) = / (-lnu)°"^ du, 



for better spacing. Many of the evaluations in |2j follow this form. Section 4.215 
in [5] consists of four examples: the first one, 4.215.1 is (|4.ip itself. The second 
one, labeled 4.215.2 and written as 

dx TT 

{-inxf ^ rv)' 

is evaluated as r(l — /i) by (|4.ip . The identity 



(4-3) / ^ = "FTTTCOsec /^tt, 



(4.4) n^^ni-f,)^-. — 

smTT/x 

yields the given form. The reader will find in [T] a proof of this identity. The section 
concludes with the special values 



/•I z:^ 

/ \/— In X dx — —— , 
Jo 2 



(4.5) 

as 4.215.3 and 4.215.4 

(4.6) 

Jo V^lna; 

Both of them are special cases of (|4.ip . 



1 da; 



The reader should check the evaluations 4.269.3 
(4.7) 



x^ ^ V— In a; dx = 
Jo 2 V P'^ 



(4.8) 
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and 4.269.4: 

xP-^ dx _ f¥ 
Jq V-lnx V P 
by reducing them to ((27T|) . Also 4.272.5, 4.272.6 and 4.272.7 

(4.9) Hilr^xf^ = T{l+p), 

Ji 

1 1 

i-lnx)"''^ x^-^dx = — r(M), 

can be evaluated directly in terms of the gamma function. 
Differentiating (|4.ip with respect to a yields 4.229.4 in [2]: 

(4.10) / ln(-lna;) (-lnx)°^^ dx = r'(a) = V(a)r(a), 
with ipla) defined in (|2.18p . The special case a = 1 is 4.229.1: 

(4.11) [ ln(-lna;) dx = -7, 

Jo 

and 

dx 

(4.12) / ln(-ln.x) " = -(7 + 21n2)V^, 
Jo V—lnx 

that appears as 4.229.3, is obtained by using the values F (i) = and (^) = 
-(7 + 21n2). 

The same type of arguments confirms 4.325.11 

fi rt^^^ dr fir 

(4.13) / ln(-lnx)^=^ = -(7 + ln4^)./-, 
Jo V-lnx V M 

and 4.325.12: 

/•I 1 

(4.14) / ln(-ln.x) (-Inx)''"^ x'^-ida;^ — ^(^)[?/'(A^)-lnH■ 
Jo 

In particular, when fi ^ 1 we obtain 4.325.8: 

/•I 1 

(4.15) / ln(-ln.T)a;''"ida; = (7 + lnty). 

Jo V 

5. The presence of fake parameters 

There are many formulas in [2] that contain parameters. For example, 3.461.2 
states that 

and 3.461.3 states that 

(5.2) r x^'^+^e-P-" dx 



2p 



n+l ' 



GAMMA FUNCTION 



9 



The change of variables t = px^ eliminates the fake parameter p and reduces 3.461 .2 
to 



(5.3) / t"-2e-*dt = 



1 



(2n- 1)!! 
2" 



and 3.461.3 to 

(5.4) / re"*(ii = n!. 





These are now evaluated by identifying them with r(?T. + i) and r(n + 1), respec- 
tively. 

A second way to introduce fake parameters is to shift the integral (|2.ip via 
s = t + b to produce 

/>oo 

(5.5) / (s-6)°-ie-''^ds = /i-°e-'^'T(a). 

Jb 

This appears as 3.382.2 in [2]. 

There are many more integrals in [2i that can be reduced to the gamma function. 
These will be reported in a future publication. 
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